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ABSTRACT 

This  paper  describes  a  multi-resolution,  distributed  design  and  simulation  software  for  rapid  prototyping  and  analysis  of 
complex  systems  using  a  Co-simulation  approach.  The  current  focus  of  this  work  is  on  the  modeling  of  the  engine  cooling  system 
in  the  Ford  Escape  Hybrid  SUV  vehicle.  In  that  particular  vehicle  the  cooling  system  consists  of  three  coupled  subsystems:  a) 
engine  cooling;  b)  electronic  transaxle  cooling;  and  battery  cooling.  This  paper  discusses  two  aspects  of  this  work:  a)  high  level 
description  of  the  developed  models  and  co-simulation  approaches;  and  b)  comparison  of  co-simulation  to  test  data.  Most  of  the 
model  predictions  deviated  from  the  test  data  by  less  than  5%.  Results  indicate  that  distributed  multi-resolution  simulations  can 
significantly  accelerate  the  analysis  of  flow-thermal  processes  in  complex  vehicle  systems.  Moreover,  the  approach  allows 
coupling  of  different  codes  with  different  functionalities  to  obtain  integrated  results  not  possible  with  any  one  individual  code. 


1.  INTRODUCTION 

Many  physical  systems  can  be  of  sufficient  complexity 
(both  geometric  and  physics-related)  to  make  fully  resolved 
3D  simulations  impractical  due  to  complicated  gridding 
requirements  and  potentially  very  slow  execution  times  on 
complex  grids.  Moreover,  a  single  software  package  may 
not  have  all  the  capabilities  required  for  complete  vehicle 
thermal  analysis.  Modem  army  and  commercial  vehicles 
with  hybrid  power  systems  that  contain  multiple  heat¬ 
generating  and  rejection  components  are  one  example  of 
such  systems,  as  shown  in  Fig.  1.  This  particular  system 
shows  three  independent  but  coupled  cooling  loops:  engine, 
motor-electronics  (M/E)  module  and  the  battery  pack  using 
the  A/C  system. 

The  distributed  multi-resolution  approach  sidesteps  these 
difficulties  by:  a)  partitioning  a  complex  system  into 
interacting  components  that  can  be  represented  by  reduced 
models  of  varying  levels  of  fidelity;  and  b)  using  several 
codes  in  coupled  parallel  or  parallel/series  execution,  each 
performing  a  set  of  specific  computational  tasks  and 
exchanging  information  in  real  time.  For  generality, 
information  exchange  takes  place  with  the  aid  of  a 
simulation  environment  that  allows  inclusion  of  additional 
system  component  models  and  legacy  codes  with  minimal 
code  modifications.  Such  a  computing  approach  is  also 
known  as  grid,  or  co-simulation. 


Figure  1:  Multi-component  Cooling  System  in  a  Ford 
Escape  Hybrid 


The  examples  discussed  in  this  paper  show  the 
advantages  of  this  approach  for  analysis  of  complex 
systems.  The  multi-resolution  co-simulation  methods  allow 
the  use  of  the  time-intensive,  high-fidelity  models  only 
where  needed,  with  the  rest  of  the  system  being  modeled 
using  much  faster,  lower  resolution  approaches. 

2.  MULTI-RESOLUTION  ANALYSIS 

Multi-resolution  analysis  partitions  a  complex  continuous 
system  into  components  where  each  component  is 
represented  by  a  separate  model.  In  our  approach  we  have 
developed  a  general  wrapper  routine  that  encapsulates  a 
large  class  of  such  models  as  stand-alone  codes,  with 
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eommunieation  interfaces  that  allow  the  eodes  to  exehange 
data  with  one  another  between  suitably  defined  boundary 
and/or  volume  eonditions.  This  approaeh  allows  simulation 
of  a  large  elass  of  dynamie  systems  to  different  levels  of 
aecuracy  since  it  allows  eoupling  of  models  defined  on  OD- 
3D  domains.  Data  exehange  is  brokered  by  the  Open 
Arehitecture,  eo-simulation  environment  CoSim  that 
eonneets  the  system  components.  Fig.  2  shows  the 
eonneetion  seheme  for  a  stand-alone  A/C  system. 


Figure  2:  CoSim-Brokered  A/C  Connection  Scheme 


The  component  models  are  eonneeted  to  a  3D  CFD 
model  that  caleulates  the  under-hood  ambient  conditions 
that  regulate  heat  exehange  between  the  system  and  the 
environment.  As  indicated,  CoSim  also  allows  visualization 
of  a  component  response  history  during  run  time.  The  main 
advantage  of  conneeting  the  system  components  through  an 
independent  environment  rather  than  conneeting  them 
directly  is  that  additional  components  ean  be  added  to  the 
system  without  affeeting  the  original  conneetivity  or  data 
transfer  synchronization. 

Although  system  partitioning  is  by  necessity  problem 
specifie,  the  possibility  of  using  separate  codes  to  simulate 
the  dynamies  of  different  components  gives  this  form  of 
multi-resolution  approaeh  a  large  degree  of  latitude  in  how 
the  systems  are  partitioned.  Multi-resolution  analysis  can  be 
earried  out  on  different  levels  and  with  different 
requirements  that  may  be,  for  example: 

\. Fully-resolved:  full  3D  analysis  performed  on  the 
vehicle-package  system  using  standard  CFD. 

2.  Mixed-resolution:  eombination  of  full  2D  or  3D 
detailed/meshed  model  and  reduced-order 
representations  of  components  within  the  system. 

3 .  System-level:  assembly  of  redueed-order  models  of 
eomponents  for  rapid  generation  of  results  for  the 
entire  vehicle-package  system. 

A.  Interface  with  legacy  codes:  eaeh  code  performs  a 
different  analysis  task.  For  example,  one  eode  ean 
generate  vehicle  surfaee  heat  flux  data  (due  to  internal 
heating)  that  another  eode  will  use  to  eompute  vehiele 
surfaee  cooling  and  thermal  signature. 


5.  Selective  focusing:  using  reduced-order  models  (with 
inereased  resolution)  on  speeifie  eomponents/ 
assemblies,  while  using  regular  reduced-order  models 
for  the  remaining  eomponents.  This  feature  allows  fast 
analyses  of  parametrie  changes  in  the  ‘foeused’ 
portions,  and  is  extremely  useful  for  rapid  thermal 
prototyping,  signature  management,  and  optimization. 

2,1  Simulation  Environment  and  Data  Exchange 
The  basic  function  of  the  simulation  environment  is  to 
sehedule  data  transfer  between  the  different  eodes  that 
eomprise  the  system  of  interest.  The  eodes  do  not 
eommunieate  directly,  but  only  through  the  environment. 
The  environment  also  eontains  transformer  funetions  that 
ean  be  used  for  data  sealing,  resealing,  and  modeling  of 
eomponents,  and  contains  a  user-expandable  library  of 
eomponent  models  of  different  fidelities,  eommon  to 
vehiele  eooling  systems.  The  environment  strueture  allows 
an  arbitrary  number  of  eodes  to  be  eonneeted  in  an  arbitrary 
fashion,  and  exehange  an  arbitrary  amount  and  type  of  data. 
Data  exehange  is  executed  using  the  CORBA 
eommunieations  protocol  [e.g.  Siegel,  2007]  allowing 
platform-to-platform  data  exehange  over  a  computer 
network.  The  protoeols  are  eneapsulated  in  a  set  of 
Application  Programming  Interfaces  (APIs).  APIs  are 
eodelets  designed  to  be  appended  to  any  stand-alone 
eomponent  code,  with  minimum  modifieation  to  the  host 
eodes.  Corresponding  APIs  exist  in  CoSim.  Fig.  3  shows  a 
sehematie  of  the  data  exchange  process. 


Server 

{CoSim) 


\ 


Data 


APIs 

\ 


Client 

(Component) 


Figure  3:  Server-Client  eonnectivity  Between  the 
CoSim  Environment  and  System  Component  Models 


The  environment  exeeutes  the  following  functions: 

1.  Data  reception  and  reeeption  scheduling  (using  APIs); 

2.  Data  distribution  and  sequencing  to  proper  models; 

3.  Data  pre-proeessing  (mainly  appropriate  sealing 
operations); 

4.  Data  proeessing  (if  library  models  are  used  in 
simulations); 

5.  Data  post-processing  (appropriate  output  resealing); 

6.  Data  transmission  and  seheduling  (using  APIs). 


Data  is  exehanged  at  a  rate  speeified  by  the  user.  Any 
data  can  be  exehanged  (pressure,  temperature,  etc.)  from 
any  point  in  the  eomputational  domain.  The  user  is 
responsible  for  providing  the  software  front-ends  that  will 
extraet  (impose)  the  required  data  from  (to)  the  code  of 
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interest  and  supply  the  data  to  the  APIs  for  transmission  (or 
reeeption).  When  eomponent  models  use  implieit  solvers, 
information  is  typically  passed  every  several  iterations  to 
make  the  simulation  tightly  eoupled.  During  the  multi- 
eomponent,  multi-resolution  simulation,  the  environment 
ensures  that  the  exeeutions  of  each  code  and  data  exehange 
are  properly  synchronized.  The  aetual  data  exehange 
proeess  is  illustrated  in  Fig.  4  for  the  eoupling  of  an 
example  radiator  model  and  the  environment  in  whieh  it  is 
embedded. 


This  approach  is  used  for  eommunieations  of  any  reduced 
model  with  the  assoeiated  CFD  domain.  A  placeholder 
domain  is  defined  in  the  CFD  grid  with  whieh  the  reduced 
model  is  exehanging  information:  the  CFD  domain  sends 
the  ambient  conditions  data  to  the  redueed  model  that  in 
turn  uses  this  (and  other)  information  to  calculate  the  (e.g.) 
heat  flux  that  is  sent  baek  to  the  CFD  solver. 

2.2  Application  to  Cooling  of  Ford  Escape  Subsystems 

Fig.  1  shows  the  three  major  sources  of  heat,  eaeh  of 
whieh  with  its  own  eooling  loop:  gasoline  engine,  M/E  unit, 
and  the  battery  paek.  It  was  possible  to  aecurately  measure 
heat  rejeetion  by  the  various  components  in  the  eooling 
systems  in  each  separate  loop,  using  similar  procedures.  For 
eaeh  loop,  working  fluid  (eoolant  or  refrigerant)  mass  flow 
and  temperature  change  aeross  eaeh  component  was 
measured  using  turbine  flow  meters  and  thermocouples, 
respeetively.  Air  velocities  through  the  heat  exehangers 
were  measured  directly  and  mass  flow  rates  were  estimated 
from  temperature  measurements  using  arrays  of 
thermoeouples.  Radiator  fan  speeds  were  measured 
optically.  The  following  general  types  of  measurements 
were  performed. 

Engine  Cooling:  was  modeled  by  eonneeting  the 
following  eomponents  together:  M/E  unit  (transmission); 
pump;  radiator,  piping,  as  well  as  transmissionEnvironment 
(for  setting  load/speed  parameters)  and 
radiatorEnvironment  (for  speeifying  data  exehange  with  the 
CFD  underHood  model).  The  measurements  ineluded: 


coolant  mass  flow,  temperature  measurements  across  the 
engine,  temperature  and  air  mass  flow  across  the  radiator. 


Figure  5:  Modeling  the  Engine  Cooling  Foop 


Fig.  5  shows  the  configuration  of  this  system,  along  with 
the  information  flow.  We  have  performed  a  number  of 
system-level  demonstration  simulations  that  show  the 
eomponent  and  system  response.  Fig.  6  shows  an  example 
of  transient  response  to  a  perturbation  to  vehiele  speed. 


It  ean  be  seen  that  all  caleulations  are  stable,  and  the 
pump  is  operating  near  its  upper  limit:  inereases  in  the 
vehiele  speed  have  little  effeet  on  the  resultant  pumping 
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rate.  The  opposite  (not  shown)  is  true  as  the  vehiele  speed 
deereases:  eoolant  flow  rate  decreases  notably,  as  expeeted. 

M/E  Unit  Cooling:  We  have  again  conneeted  six 
eomponents  that  form  the  transaxle  eooling  subsystem.  The 
eooling  eonfiguration  is  shown  in  Fig.  7. 


Figure  7:  Modeling  the  M/E  Unit  Cooling  Eoop 

The  eomponents  and  the  information  flow  are  similar  to 
that  for  the  engine,  replacing  the  latter  with  transaxle 
model,  the  engine  pump  with  transaxle  pump,  and  the 
engine  radiator  eonfiguration  with  one  appropriate  for 
transaxle  eooling.  The  measurements  were  the  same  as  for 
the  engine:  coolant  mass  flow,  temperature  measurements 
across  the  transaxle,  temperature  and  air  mass  flow  across 
the  radiator.  This  system  is  also  strongly  coupled  sinee 
ehanges  in  the  radiator  outlet  eoolant  temperatures  directly 
affect  the  temperature  inerease  aeross  the  transmission. 

As  for  the  engine,  we  have  performed  a  number  of 
system-level  demonstration  simulations  that  show  the 
eomponent  and  system  response.  Fig.  8  shows  an  example 
of  transient  response  to  a  perturbation  to  vehiele  speed.  We 
see  the  same  mass  flow  trends  that  we  have  observed  for 
load  variation:  notieeable  eoolant  flow  increase  with 
inereasing  vehicle  veloeity,  and  smaller  flow  increase  with 
deereasing  velocity.  Again,  calculations  are  stable.  Note  the 
slight  increase  in  exit  air  veloeity  due  to  air  heat-up,  as 
shown  in  panel  (e).  We  note  that  although  the  radiator 
dimensions  for  the  transaxle  and  engine  radiators  are 
similar,  ICE  engine  eoolant  flow  rate  is  more  than  an  order 
of  magnitude  higher  than  in  the  transaxle  eooling 
subsystem,  sinee  the  former  generates  considerably  more 
heat  than  the  latter.  This  inereased  heat  results  in 
eonsiderably  higher  temperatures  and  radiator  powers,  than 
for  the  transaxle.  System  responses  to  variations  in  the 
vehicle  speed  conditions  are  also  much  less  pronounced. 


Battery  Pack  Cooling:  Fig.  9  shows  the  eonfiguration  of 
this  system,  along  with  information  flow. 


The  system  was  modeled  by  conneeting  the  following 
eomponents  together:  eondenser,  evaporator,  compressor, 
expansion  valve,  piping,  as  well  as  the  vehieleEnvironment 
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(for  setting  A/C  operating  parameters)  and  evaporator 
CondenserEnvironment  (for  speeifying  data  exehange  with 
the  CFD  underHood  model). 

As  for  the  engine  and  the  transaxle,  we  have  performed  a 
number  of  system-level  demonstration  simulations  that 
show  the  component  and  system  response.  Fig.  9  shows  an 
example  of  transient  evaporator  pressure  and  temperature 
responses  to  perturbations  in  the  refrigerant  flow  rates. 
Once  again,  all  the  caleulations  are  stable. 


2.3  Application  to  Cooling  of  Ford  Escape  System 
The  eomplete  eooling  system  consists  of  the  three 
subsystems  diseussed  previously.  To  demonstrate  stable 
interactions  with  a  more  realistie  environment,  we  have 
eonneeted  the  four  heat  exehangers  (two  radiators, 
eondenser  and  the  evaporator)  to  the  STAR-CCM+  CFD 
eode  that  simulates  the  3D  thermal-flow  fields  in  very 
simplified  underHood  (for  heat  rejeetion)  and  battery 
eompartment  (for  heat  absorption)  domains.  Fig.  10  shows 
the  complete  system  as  represented  by  the  interconneeted 
subsystems.  We  note  that  the  following  results  are 


qualitative  in  nature  since  the  simulations  were  performed 
using  a  very  simplified  eomputational  domain. 


ME  (TransAxle)  Subsystem 


Engine  Subsystem 


CFD  Environment 


UnderHood 

Compartment  Battery 

(’ompartment 


t 


C(m4tn%fr 


Environment-starCCM  CFD 


hvaporafttr 


AC  Subsystem 


(  Cooling  —  Heating ) 

Figure  10:  Sehematie  of  the  eomplete  system 


The  CFD  domain  was  divided  into  two  isolated 
subDomains,  underHood  and  the  battery  compartment,  eaeh 
using  different  sources  of  external  air.  The  two  do  however 
eommunieate  indirectly  through  the  refrigerant  in  the  A/C 
system  that  heats  up  through  battery  eooling,  and  eools 
down  through  heat  rejection  in  the  underHood  domain.  Two 
mesh  models  and  two  physieal  models  were  seleeted  for 
those  two  domains,  respeetively.  Unstructured  hexahedron 
grids  were  generated  using  the  3D-CAD  faeility  in 
starCCM+. 

Fig.  11a  shows  the  eomposite  environment  indieating  the 
embedded  heat  exehangers  and  the  respective  air  flows.  In 
this  eo-simulation,  starCCM+  provided  the  air  velocity, 
temperature  and  pressure  to  the  four  heat  exchangers.  The 
direction  of  the  air  velocity  vectors  determined  the  direetion 
of  data  exchange:  the  front  faee  was  in  the  upwind  direetion 
and  the  rear  faee  was  downwind.  CFD  data  on  the  front  faee 
of  the  heat  exchangers  was  mapped  onto  the  local 
exchanger  face  boundary,  and  was  then  passed  to  the 
exehanger  models.  The  models  then  caleulated  a  new 
thermal  state  of  the  exiting  air  at  the  rear  faee,  whieh  was 
then  mapped  onto  the  corresponding  position  in  the  CFD 
domain.  The  grid  mapping  proeedures  were  used  on  the  8 
interacting  faces  (4  eomponents  times  2  faees/eomponent) 
that  eonneeted  the  CFD  and  the  heat  exchanger  model 
domains. 
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Unsteady  simulations  were  performed  in  which  data 
exchange  between  the  system  components  was  executed  at 
every  time  step,  with  100  iterations  per  time  step.  Fig.  lib 
shows  temperature  distributions  at  a  seleeted  timestamp,  in 
the  two  eompartments  in  the  middle  cutting-plane.  The  case 
was  based  on  the  operations  of  the  Ford  Eseape  test  vehiele, 
using  the  Escape’s  eomponent  geometries  with  the 
following  compressor  and  valve  eoeffieients:  compressor 
speed  =  3000  RPM;  valve  area  control  parameter  =  17.6. 
We  note  that  the  underHood  and  the  battery  eompartments 
are  grossly  simplified  sinee  the  intent  of  the  simulations  is 
not  to  predict  the  cooling  characteristics  of  this  particular 
test  vehiele,  but  to  show  the  eapabilities  for  mixed- 
resolution  simulations  of  vehiele  (or  other)  eooling 
systems.  It  ean  be  seen  that  the  underHood  heat  exehangers 
(transaxle  radiator,  eondenser  and  ICE  engine  radiator)  heat 
the  environment,  whereas  the  battery  heat  exchanger 
(evaporator)  eools  the  battery  domain.  The  results  are  as 
expeeted  and  are  physically  reasonable.  For  this 
eonfiguration,  there  is  little  differenee  between  steady  and 
unsteady  simulations. 


3.  Data  Comparison 

This  seetion  doeuments  the  performance  of  the  redueed 
models  with  eomparisons  to  vehiele  test  data.  The  data  were 
generated  from  35  measurements  made  at  7  different  speeds 
(10  mph,  20  mph,  30  mph,  40  mph,  50  mph,  60  mph  and  70 
mph)  and  5  different  vehiele  loads.  Eaeh  run  lasted  60 
seeonds  (except  for  30  seconds  at  wide  open  throttle)  and  a 
measurement  was  made  every  seeond.  For  eaeh 
measurement,  120  different  physieal  quantities  were 
aequired.  Therefore  7  X  5  X  120  x  60  =  252,000  raw  data 
points  were  recorded.  As  will  be  seen  in  the  following, 
eonsiderable  post-proeessing  was  required  in  order  to  make 
meaningful  eomparisons  between  the  models  and 
experiments.  In  some  cases  data  scatter  was  signifieantly 
high,  greater  than  the  average  values  of  the  measured 
quantities.  In  other  eases,  eomparison  data  had  to  be 
inferred.  For  this  inferred  data,  the  following  two  variables 
assoeiated  with  the  A/C  system  were  of  main  concern. 

3.1.1  Temperature:  evaporator.  These  data  were  taken 
only  aeross  the  expansion  valve,  but  not  at  the  evaporator 
outlet.  This  quantity  was  inferred  from  other  measurements 
as  follows.  We  estimated  the  inlet  enthalpy  from  the  inlet 
valve  temperature,  estimating  the  pressure  and  noting  that: 
a)  the  liquid  enthalpy  varies  very  little  with  pressure,  and  b) 
flow  aeross  the  valve  is  essentially  isenthalpic.  The  quality 
was  estimated  from  the  inlet  temperature  and  enthalpy, 
using  the  lever  rule  for  enthalpy  in  the  two-phase  region 
(that  is  =  xh^+{\- x)h^\  x  =  fluid  quality).  We 

assumed  that  the  inlet  enthalpy  to  the  evaporator  is  equal  to 
the  outlet  enthalpy  to  the  valve.  The  evaporator  outlet 
temperature  was  estimated  by  noting  that  the  valves  are 
typically  tuned  so  that  the  temperature  differenee  across 
evaporators  is  in  the  range  5-10  °C.  Furthermore,  the  bulk 
of  the  enthalpy  received  by  the  evaporator  oecurs  in  the 
two-phase  region  and  thus  the  enthalpy  received  in  the 
gaseous  region  is  small  eompared  to  the  former.  Thus,  the 
error  in  estimating  the  evaporator  outlet  enthalpy  remains 
relatively  small.  For  all  evaporator  eomparisons,  we  ehose 
the  midway  point  so  that  the  refrigerant  at  the  outlet  is 
heated  by  7  °C  over  the  inlet.  Data  aeross  the  condenser 
were  well  defined. 

3.1.2  Pressure:  both  the  evaporator  and  condenser. 

1  Evaporator  pressure  was  estimated  from  the  inlet 
refrigerant  temperature,  noting  that  pressure  depends 
only  on  temperature  that  does  not  ehange  appreeiably. 

2  Condenser  pressure  was  estimated  at  1.6  MPa,  based 
on  typical  condenser  operations.  This  pressure  is 
required  for  speeification  of  the  inlet  enthalpy  of  the 
refrigerant  that  is  in  the  gas  state  and  varies  with  both 
temperature  and  pressure. 
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3.1.3  Data  Consistency  Checks —  a  more  serious  problem 
with  experimental  comparisons  was  related  to  the  quality  of 
test  data,  and  in  partieular  to  data  for  the  heat  exehangers. 
Closer  examination  showed  that  the  measured  data  was 
ineonsistent,  in  the  sense  that  radiator  powers  measured 
from  airside  and  liquid-side  (eoolant)  measurements  did  not 
mateh  within  a  reasonable  range.  Here,  heat  exehanger 

power  P  is  defined  by  E  „  =  m  „  C  5T  ,  . 

^  cool  fair  cool!  air  p, cool  I  air  cool  fair 

Explieitly,  eoolant  and  airside  powers  are  P  =  m  C^  ST  , 

P  =  p  Av  C  ST  .  The  differenees  (dT  ,  5T  )  were 

measured  at  front  and  baek  faees,  and  coolant  inlets  and 
outlets,  respeetively.  Under  ideal  conditions  we  expect  that 

P=P.  However,  many  measurement  eonditions  were  far 
from  ideal  and  very  often  measurements  were  reeorded  for 
whieh  Pi^P.  Thus  as  an  aid  to  data  interpretation,  all  the 

eomparisons  for  heat  exchanger  data  inelude  the  power 
eonsisteney  ehecks. 

3.1.4  Averaging  of  Values  and  Other  Uncertainties — 
although  both  the  experiments  and  models  use  spatially 
loeal  data,  final  eomparisons  were  based  on  spatially 
averaged  data  for  the  heat  exehanger  airside  temperatures 
and  air  veloeities.  Tests  however  show  high  non¬ 
uniformities  in  both  velocity  and  temperature  fields  at  inlet 
and  outlet  faees  of  the  heat  exehangers.  Although  such 
nonuniformities  are  very  likely  to  affeet  the  performanee  of 
the  devices,  we  did  not  take  them  into  aecount.  For 
simplieity  we  used  an  unweighed  averaging  proeess 
although  tests  showed  that  some  data  was  signifieantly 
skewed/biased.  Two  additional  uncertainties  are  expeeted: 

1  Flow  Data:  Under  some  operating  eonditions, 
refrigerant  flows  were  in  the  laminar-turbulent 
transitions  region.  For  sueh  regions,  heat  transfer 
eoeffieients  (erucial  for  model  performanee)  are  not 
very  reliable. 

2  Heat  Exchanger  Geometry:  Although  the  overall 
external  dimensions  (e.g.  eore  height,  width,  depth, 
number  of  tubes)  were  known,  their  internal  structure 
was  not.  In  particular  the  following  had  to  be  assumed: 
Tube  arrangement  (parallel,  or  serpentine),  chose 
parallel,  with  1  tube/stride;  Louvers:  ehose  10 
louvers/fm;  inelined  at  28°,  height  =  0.1mm;  Internal 
microchannels  were  speeified  to  be  the  same  as  the 
number  of  louvers;  Tube  thickness:  speeified  E5  mm; 
Internal  wall  thickness:  speeified  0.3mm. 

Despite  the  above  uneertainties,  model  predietions  were 
remarkably  aceurate  in  comparison  to  test  data,  and  most 
fell  within  the  required  5%  difference. 


3.2  Comparison  of  Generic  Models 

Generie  models  refer  to  formulations  with  no  speeifie 
referenee  to  any  partieular  vehiele.  In  this  ease  they  refer  to 
all  single  and  multi-phase  finned  heat  exehanger  models. 

3.2.1  Engine  Radiator  Data  Consistency —  The  engine 
radiator  geometry  is  given  by:  Core  Area:  698.5  mm  x 
469.9  mm;  Tube  Geometry:  15.6  mm  x  2.22  mm;  Tube 
Pitch:  41  tubes  @17  mm/tube.  Fig.  12  compares  the  air 
power  Pa  with  the  coolant  power  Pc  for  the  engine  radiator. 

We  show  {P  -P  )IP^  expressed  in  pereentages. 
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Figure  12:  Relative  differences  between  measured 
air  and  engine  coolant  powers 

The  measured  air-based  power  estimates  are 
approximately  30%  higher  in  comparison  to  coolant-based 
estimates.  One  of  the  reasons  for  this  discrepancy  may  be 
due  to  the  positioning  of  the  air-velocity  probe  loeated 
between  the  eondenser  and  the  engine  radiator.  For  probe 
loeations  sufficiently  far  from  the  radiator  face,  air  veloeity 
going  through  the  radiator  fins  would  be  overestimated  due 
to  partial  bloekage  of  the  eore  by  the  internal  radiator 
strueture  (fins,  tubes  and  louvers).  This  bloekage  would 
have  slowed  down  the  air  going  through  the  radiator,  and 
eaused  some  air  to  spill  out  the  sides  of  the  radiator.  Sinee 
we  did  not  correet  the  airflow  readings  for  such  blockage 
effects,  airspeed  inside  the  radiator  is  likely  to  be 
overestimated,  aceounting  for  the  higher  air  power  values. 

3.2.2  Engine  Radiator  Model  Predictions —  Fig.  13 
eompares  the  relative  differences  in  coolant  power  between 
model  and  experiments,  for  a  41 -tube  radiator.  The 
microehannel  area  tuning- faetor  was  set  to  1 .  Best  fit  to 

the  experiment  is  obtained  for  tuning  factor  -  1.7  .  We 

note  that  in  the  laminar  region  (small  mass  flow  rate),  this 
difference  is  in  the  range  of  +10%  with  little  seatter,  while 
in  the  turbulent  region,  the  differenee  is  in  the  range  of  - 
10%  with  large  seatter.  These  differenees  in  the  model 
prediction  can  be  explained  by  the  faet  that  we  use  two 
different  heat  transfer  eorrelations  for  laminar  and  turbulent 
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flows  and  by  the  fact  that  the  measurements  are  more 
impreeise  for  high  coolant  mass  flow  (see  below). 


Heat  transfer  x  1 .7 

Figure  13:  Comparison  model/experiment  data  for 
engine  coolant  power  for  heat  transfer  faetor=1.7 

We  note  the  following:  during  the  measurements,  air 
velocity  has  been  essentially  constant  at  approximately  1 
m/s,  while  the  coolant  mass  flow  varied  from  a  run  to 
another  by  2  orders  of  magnitude  between  30  g/s  and  2  kg/s. 
Therefore  when  the  coolant  mass  flow  is  high,  air  takes  out 
little  energy  per  amount  of  coolant  and  therefore  the 

difference  5T  =T  -  T  is  small.  Measurements  were  taken 

c  2  1 

of  only  and  but  not  the  temperature  difference. 

Therefore,  the  relative  uneertainty  on  the  measurement  of 
the  eoolant  temperature  differenee  will  be  large  (on  the 
order  of  1 00%  in  some  measurements)  when  the  uneertainty 
in  measured  temperatures  is  large.  This  fact  explains  why 
the  measured  seatter  on  the  relative  differenee  between 
model  and  experiment  is  larger  for  high  mass  flows.  Also 
for  high  mass  flows  coolant  flow  is  turbulent  while  for  low 
mass  flows  it  is  laminar.  The  eorrelations  used  to  compute 
the  Nusselt  number  Nu  differ  significantly  in  these  two 
regimes.  For  laminar  flow,  Nu  is  a  eonstant  that  only 
depends  on  the  geometry  of  the  tube  (approximately  2.9) 
while  for  turbulent  mass  flows,  the  Nusselt  number  depends 
on  physical,  geometrical  and  flow  parameters  (see 
Gnielinski  [5])  Between  the  two  flow  regimes  there  exists 
an  ill-defined  transition  zone  that  oecurs  at  approximately 
500  g/s.  This  area  has  been  represented  in  green  in  the 
figures  above,  and  it  is  in  this  area  that  most  of  the  model 
predictions  deviate  from  the  test  data. 

3.2.3  Transaxle  Radiator  Data  Consistency —  The 
transaxle  radiator  geometry  is  given  by:  Core  Area:  730.25 
mm  X  304.8  mm;  Tube  Geometry:  16.13  mm  x  2.1  mm; 
Tube  Pitch:  26  tubes  @  28  mm/tube.  Fig.  14  eompares  the 
air  power  Pa  with  the  eoolant  power  Pc  for  the  engine 
radiator. 


Figure  14:  Relative  differences  between  measured 
air  and  transaxle  coolant  powers 


The  figure  shows  that  the  data  appears  to  be  distributed 
into  three  groups.  We  do  not  have  a  elear  explanation  for 
this  pattern:  In  the  first  set,  air  power  is  underestimated  by 
25%;  in  the  second,  air  power  is  underestimated  by  only 
5%;  the  third  corresponds  to  the  outliers  mentioned  below. 

3.2.4  Transaxle  Radiator  Model  Predictions —  The 
transaxle  radiator  analysis  (Fig.  15)  is  similar  to  that 
performed  for  the  engine  radiator.  We  note  that  in  all  tests 
and  simulations,  the  coolant  mass  flows  are  almost  eonstant 
and  are  considerably  smaller  than  those  in  the  engine 
radiator:  they  are  in  the  range  of  only  135  g/s  to  150  g/s. 


Figure  15:  Comparison  model/experiment  data  for 
transaxle  eoolant  power,  heat  transfer  faetor=  2.6 

These  small  mass  flows  result  in  laminar  flow  conditions 
for  which  the  Nusselt  number  is  constant  at  approximately 

2.9.  The  best  fit  to  experiment  is  obtained  for  =  2.6  .  For 

both  cases,  most  of  the  eomparisons  show  that  the  model 
matches  the  experiment  within  the  required  5%.  Also, 
several  model  predietions  at  higher  mass  flows  overestimate 
the  power  by  20%  (or  underestimate  it  by  10%).  This  trend 
seems  to  be  assoeiated  with  the  outliers  shown  in  Fig.  15.  It 
is  possible  that  during  the  measurements,  the  coolant 
temperature  differenee  may  have  been  overestimated  for 
these  points. 
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3.2.5  Evaporator  Data  Consistency —  The  evaporator 
geometry  is  given  by:  Core  Area:  203.2  mm  x  127  mm; 
Tube  Geometry:  57.15  mm  x  1.7  mm;  Tube  Pitch:  10  tubes 
@  12.7  mm/tube.  Fig.  16  compares  air  power  Pa  to 
refrigerant  power  Pr  for  the  evaporator  for  full  and  reduced 
test  data  sets,  respectively. 


_ ^ ^ _ I _ I _ ^ ^ _ I _ ^ ^ _ 
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Figure  16:  Comparison  air/refrigerant  measured 
powers  (evaporator  eomplete  raw  data) 

As  diseussed  at  the  beginning  of  this  Seetion  3.1,  there 
exist  problems  with  the  integrity  of  data  associated  with 
A/C  system  operations.  In  the  comparisons  we  define  air- 

based  power  Pa  ealculations  using  P  =  p  Av  C  ST  .  The 

^  ^  a  ’  a  a  pa  a 

difference  (5r  is  well  defined  sinee  air  temperatures  were 
measured  in  front  and  behind  the  evaporator.  The  power 
P=m{h^^-h,J  is  not  well  defined  since  refrigerant 

temperature  has  been  measured  across  the  expansion  valve, 
and  not  across  the  evaporator.  However,  based  on  the 
discussion  at  the  beginning  of  this  seetion,  these  values  ean 
be  estimated  with  a  reasonable  aceuraey  based  on  well- 
chosen  assumptions.  We  ean  estimate  the  inlet  enthalpy  hin 
and  derive  its  saturated  pressure  Psat  from  the  fact  that  the 
fluid  is  in  two-phase  state.  This  is  the  refrigerant  pressure 
inside  the  evaporator. 

Fig.  16  shows  power  eomparisons  for  all  the  measured  35 
test  points.  It  ean  be  seen  that  for  most  of  the  data  the 
measured  air-based  power  estimates  are  eonsiderably 
smaller  than  the  estimated  refrigerant-based  powers  that 
depend  on  the  refrigerant  mass  flow  rate.  In  some  cases,  air 
powers  are  two  orders  of  magnitude  smaller  than  those  of 
the  refrigerant.  In  such  cases,  the  relative  differenee 
between  air  and  refrigerant  powers  is  close  to  -100%:  a 
loeus  represented  by  a  red  line  in  the  figure.  This  means  that 
in  praetice,  the  evaporator  is  not  eooling  the  air.  The  mass 
flow  rate  is  derived  from  the  volumetric  flow  rate  by 
assuming  that  the  refrigerant  is  liquid.  We  conelude  that  in 
these  cases,  the  volumetric  flow  rate  was  measured  when 


the  refrigerant  was  gaseous,  and  therefore,  the  air 
conditioning  system  was  not  working  properly. 

In  order  to  obtain  data  suitable  for  model-test 
eomparisons,  we  pruned  the  test  data  by  defining  a  cut-off 
ratio  P^Ip^=1.5,  beyond  whieh  the  A/C  system  was 

deemed  inoperative.  This  pruning  redueed  the  number  of 
valid  test  data  points  to  10. 
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Figure  1 7:  Comparison  air/refrigerant  measured 
powers  (evaporator,  pruned  valid  data) 

The  effects  of  this  pruning  are  shown  in  Fig.  17  that 
eompares  the  air  power  with  the  refrigerant  power  for  the 
valid  cases.  Still,  the  remaining  data  is  not  perfect  since  at 
low  mass  flow  rates  the  air  power  is  overestimated,  while  at 
high  mass  flow  rates  the  air  power  is  underestimated.  This 
discrepancy  is  likely  due  to  the  thermal  inertia  of  the 
evaporator. 
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3.2.6  Evaporator  Model  Predictions —  Fig.  18  compares 
the  relative  differences  in  eoolant  power  between  model  and 
experiment  for  an  evaporator. 
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Heat  transfer  faetor  =  2 

Figure  18:  Comparison  model/experiment  for  evaporator 
refrigerant  power  (fin  efficieney  =  0.43) 


In  all  cases  the  mieroehannel  tuning  eoeffieient  was 
set  to  1,  with  the  insulated  fin  approximation,  or  efficieney 
=  0.43.  Best  fit  to  data  is  obtained  for  2.  For  this 
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value  most  of  the  data  lies  within  2%  and  all  data  lies 
within  5%.  Sueh  good  aceuracy  between  the  model  and  the 
experiment  ean  be  explained  by  the  faet  that  most  of  the 
energy  transfer  between  the  air  and  the  refrigerant  oecurs 
when  the  fluid  is  evaporating  and  requires  large  amounts  of 
latent  heat.  In  the  first  approximation,  the  energy  flux 

P  ~  rhx  h^^  .  The  amount  of  energy  transferred  in  the  gas 

region  is  comparatively  small  and  therefore  the 
measurement  error  of  the  inlet  and  outlet  temperatures  will 
have  little  impaet  on  the  measured  enthalpy  differenee. 


Fig.  20  that  compares  the  air  power  with  the  refrigerant 
power  for  the  valid  cases  using  pruned  data.  Again,  the 
remaining  data  is  not  perfeet  sinee  at  lower  mass  flow  rates 
the  air  power  is  overestimated,  while  at  higher  mass  flow 
rates  the  air  power  is  underestimated.  As  for  the  evaporator, 
this  diserepaney  is  likely  due  to  the  thermal  inertia  of  the 
eondenser. 

3.2.8  Condenser  Model  Prediction  Comparisons —  Fig.  21 
eompares  the  model-test  refrigerant  powers  for  the  1 0  eases 
where  the  A/C  was  working. 


3.2. 7  Condenser  Data  Consistency —  Condenser  geometry 
is  given  by:  Core  Area:  654.  mm  x  419.1  mm.  Tube 
Geometry:  20.32  mm  x  2.18  mm;  Tube  Pitch:  41  tubes  @ 
10.2  mm/tube.  As  for  the  evaporator.  Fig.  19  eompares  air 
and  Pa  refrigerant  power  Pr  for  the  condenser  for  full  and 
redueed  test  data  sets. 
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Figure  19:  Comparison  air/refrigerant  measured 
powers  (condenser,  complete  raw  data) 

It  is  seen  that  for  the  full  data  set  (and  as  for  the 
evaporator)  air-based  powers  are  overestimated  for  low 
mass  flow  rates.  This  is  likely  due  to  the  faet  that  the  body 
of  the  eondenser  does  not  eool  down  instantaneously  and 
even  if  the  A/C  system  stops  working,  the  external  air  will 
still  cool  down  the  condenser. 


Figure  20:  Comparison  air/refrigerant  measured 
powers  (eondenser,  pruned  valid  data) 
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Figure  21:  Comparison  model/experiment  for  condenser 
refrigerant  power;  heat  transfer  faetor  =  0.8 

The  figure  shows  that  =  0.8  is  the  best  mateh  between 

the  test  and  model  results.  Exeept  for  the  case  where  the 
estimated  refrigerant  power  is  9500  W,  the  model  and  the 
tests  match  within  5%.  It  is  not  clear  why  our  model 
underestimates  the  refrigerant  power  for  high  volumetrie 
flow  rates.  One  possible  explanation  is  that  the  volumetric 
flow  rate  was  measured  when  the  refrigerant  was  not  in 
fully  liquid  state.  In  this  ease,  the  mass  flow  rate  would  be 
overestimated  and  therefore  the  measured  refrigerant  power 
would  be  overestimated  as  well. 

3. 3  Comparison  of  Ford  Escape  Models 

These  models  were  developed  and  tuned  for  one  speeific 
vehiele.  They  are  all  based  on  Artifieial  Neural  networks 
ANNs  trained  on  data  generated  from  the  vehicle  tests.  The 
eomparisons  show  plots  of  predieted  vs.  training  data, 
ineluding  experimental  seatter. 

3.3.1  Direct  ANN  Engine  model —  was  trained  to  directly 
correlate  the  engine  coolant  temperature  to  the  engine  state 
given  by:  vehiele  speed  (m/s),  load  (N),  coolant  volumetric 
flow  rate  (Ipm)  and  inlet  eoolant  temperature  (°C).  The 
ANN  was  speeified  in  statie  eonfiguration  with  a  4-5-1 
topology  (4  inputs,  5  neurons  in  one  hidden  layer  and  1 
output).  The  model  was  trained  such  that  the  normalized 
eumulative  predietion  error  deereased  by  2  orders  of 
magnitude. 
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Figure  22:  ANN-Engine  model  prediction  of  outlet 
temperature  (°C) 


Model-test  data  eomparisons  are  shown  in  Fig.  22.  The 
figures  show  that  the  ANN  model  predietions  are  well 
within  the  desired  5%  of  test  data  (within  1%  in  this  ease). 
3.3.2  Coolant  Pump  Models —  relate  the  coolant  mass  flow 
rates  in  the  engine  and  transaxle  cooling  loops  to  the  vehiele 
speed  (m/s)  and  load  (N).  Both  static  ANNs  were  speeified 
with  a  2-5-5- 1  topology  (2  inputs,  5  neurons  in  first  hidden 
layer,  5  neurons  in  seeond  and  1  output).  The  results  are 
shown  in  Figs.  23  and  24,  and  show  that  all  ANN  model 
predietions  are  within  the  desired  5%  of  test  data. 


Figure  23  :  ANN-Pump  model  predietion  of  engine 
coolant  flow  rate  (1pm) 


Figure  24:  ANN-Pump  model  prediction  of  transaxle 
coolant  flow  rate  (1pm) 


3.3.3  Transaxle  Heating  Model —  relates  the  coolant 
temperature  at  the  inlet  to  the  transaxle  radiator  (transaxle 
outlet)  to  the  to  the  vehiele  speed  (m/s),  load  (N),  coolant 
volumetric  flow  rate  (1pm)  and  inlet  temperature  (°C).  The 
ANN  was  specified  in  a  static  configuration  and  had  a  4-5-1 
topology  (4  inputs,  5  neurons  in  one  hidden  layer  and  1 
output).  The  model  was  trained  sueh  that  the  eumulative 
predietion  error  deereased  by  3  orders  of  magnitude. 


Figure  25:  ANN-M/E  Unit  model  predietion  of  outlet 
temperature  (°C) 


Model-test  data  eomparisons  are  shown  in  Fig.  25.  The 
figures  show  that  the  ANN  model  predietions  are  well 
within  the  desired  5%  of  test  data  (within  0.5%  in  this  ease). 

4.  SUMMARY  AND  CONCLUSIONS 

Simulation-based  prototyping  of  eomplex  systems  ean  be 
very  ehallenging  in  terms  of  domain  grid  eomplexity, 
density  and  the  assoeiated  long  simulation  times.  A 
praetieal  alternative  is  to  partition  sueh  systems  into 
interacting  component  models  of  desired  resolution  and 
fidelity  and  eouple  the  exeeution  of  these  models  into  an 
integrated  eomputational  scheme.  The  resultant  multi¬ 
resolution  eo-simulations  ean  provide  an  optimum  balance 
between  desired  aceuracy,  simulation  times,  and 
eomputational  resourees.  Moreover,  interehanging 
eomponent  models  of  different  resolution  allows  one  to 
seleetively  foeus-in  on  the  details  of  their  operation,  in  the 
eontext  of  whole-system  dynamies.  This  approaeh  is  of 
partieular  importance  in  the  virtual,  simulation-based 
prototyping  of  thermal  management  strategies  for  complex 
eooling  systems  embedded  in  vehiele  under-the-hood 
systems  for  whieh  the  use  of  full  CFD  analysis  is 
impraetical. 

The  first  two  examples  showed  the  aceuracy,  speed  and  the 
stability  of  the  multi-resolution  approaeh.  In  general,  the 
aeceleration  of  exeeution  times  eompared  to  full  CFD 
analysis  is  due  to  lower  number  of  computational  cells 
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resulting  from  replaeement  of  fully  gridded  system 
eomponents  by  equivalent  redueed  models.  Exeeution  of 
even  eomplex  redueed  models  sueh  as  those  assoeiated  with 
the  A/C  system  ean  be  timed  in  terms  of  milliseeonds;  in 
eomparison  to  their  gridded  equivalents,  the  eomputational 
time  spent  in  exeeution  of  sueh  models  is  therefore 
essentially  insignifieant.  For  modern  vehieles  with  multiple 
heat  generating  eomponents,  the  multi-resolution  approaeh 
ean  thus  be  the  enabling  methodology  for  performing 
system-level  thermal  management  prototyping  simulations. 
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